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The K+N —> K+ 5 reaction is studied for center-of-momentum energies ranging from threshold to 
3 GeV in an effective Lagrangian approach that includes the hyperon s- and u-channel contributions 
as well as a phenomenological contact amplitude. The latter accounts for the rescattering term in the 
scattering equation and possible short-range dynamics not included explicitly in the model. Existing 
data are well reproduced and three above-the-threshold resonances were found to be required to 
describe the data, namely, the A(1890), E(2030), and E(2250). For the latter resonance we have 
assumed the spin-parity of J p = 5/2“ and a mass of 2265 MeV. The E(2030) resonance is crucial in 
achieving a good reproduction of not only the measured total and differential cross sections, but also 
the recoil polarization asymmetry. More precise data are required before a more definitive statement 
can be made about the other two resonances, in particular, about the E(2250) resonance that is 
introduced to describe a small bump structure observed in the total cross section of K~ + p —> 
K + + S _ . The present analysis also reveals a peculiar behavior of the total cross section data in 
the threshold energy region in K~ + p —> K + + E _ , where the P- and D- waves dominate instead 
of the usual S- wave. Predictions for the target-recoil asymmetries of the K + N —> K + E reaction 
are also presented. 

PACS numbers: 13.75.Jz, 13.60.Rj, 13.88.+e, 14.20.Jn 


I. INTRODUCTION 

Hadron spectroscopy is an essential part of the in¬ 
vestigation to understand the non-perturbative regime 
of Quantum Chromodynamics (QCD). In principle, an 
ab-initio approach to hadron resonance physics can be 
provided by lattice QCD simulations. In particular, the 
spectra of excited baryons observed in the recent lattice 
simulations [2, 2] hold the promise of explaining the rich 
dynamics in the resonance energy region in the near fu¬ 
ture. Once quark masses drop towards more reasonable 
values and finite volume effects are fully under control, 
a close comparison to experimental data will be possi¬ 
ble. Other approaches such as the dynamical Dyson- 
Schwinger |, constituent quark models SI , and the 
Skyrme model | also generate resonance spectra. Unita- 
rized Chiral Perturbation Theory also provides a comple¬ 
mentary picture of some of the low-lying resonances |7|, i8j. 
To compare these theoretical results with the experimen¬ 
tal data, a reliable reaction theory capable of identi¬ 
fying resonances and extracting the corresponding res¬ 
onance parameters is required. Such reaction theories, 
based on a coupled-channel approach, have been devel¬ 
oped at various degrees of sophistication and are being 


improved @-111. So far, most of the experimentally ex¬ 
tracted baryon resonances come from the pion-induced 
reaction experiments, especially the nN scattering, and 
about 16 nucleon resonances and 11 A resonances have 
been identified E3- A number of A and E baryons, which 
are particles with strangeness quantum number S = — 1, 
have been also discovered E3- A review on the status of 
baryon spectroscopy is given, e.g., in Ref. jl8f 

Although the multi-strangeness baryons (S < — 1) 
have played an important role in the development of our 
understanding of strong interactions, and thus, should 
be an integral part of any baryon spectroscopy pro¬ 
gram, the current knowledge of these baryons is still ex¬ 
tremely limited. In fact, the SU(3) flavor symmetry al¬ 
lows as many S = —2 baryon resonances, called 5, as 
there are N and A resonances combined (~ 27); how¬ 
ever, until now, only eleven 5 baryons have been dis¬ 
covered E3- Among them, only three [ground state 
5(1318)1/2+, 5(1538)3/2+, and 5(1820)3/2-] have their 
quantum numbers assigned^ This situation is mainly 
due to the fact that multi-strangeness particle produc¬ 
tions have relatively low yields. For example, if there are 
no strange particles in the initial state, 5 is produced 
only indirectly and the yield is only of the order of nb 
in the photoproduction reaction [19(, whereas the yield 
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1 The parity of the ground state S has not been measured explicitly 
yet, but its assignment is based on quark models and SU(3) flavor 
symmetry. 
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is of the order of pb [20] in the hadronic A'-induced re¬ 
action, where the 2 is produced directly because of the 
presence of an S = —1 AT meson in the initial state. The 
production rates for f l baryons with S = —3 are much 
lower [[2l| . 

The study of multi-strangeness baryons has started to 
attract a renewed interest recently. Indeed, the CLAS 
Collaboration at Thomas Jefferson National Accelerator 
Facility (JLab) plans to initiate a 2 spectroscopy pro¬ 
gram using the upgraded 12-GeV machine, and measure 
exclusive f l photoproduction for the first time [22|]. Some 
data for the production of the 2 ground state, obtained 
from the 6-GeV machine, are already available jl9lL They 
were analyzed by some of the present authors [23l . |24[ 
within an effective Lagrangian approach. J-PARC is go¬ 
ing to study the 2 baryons via the K + N —> K + 2 
process (which is the reaction of choice for producing 2) 
in connection to its program proposal for obtaining in¬ 
formation on 2 hypernuclei spectroscopy. It also plans to 
study the tt + N K + K - 1- 2 reaction as well as Q pro¬ 
duction [25],[2(|. At the FAIR facility of GSI, the reaction 
p + p —y 2 + 2 will be studied by the PANDA Collab¬ 
oration [27|. Quite recently, lattice QCD calculations of 
the baryon spectra, including those of 2 and Q baryons, 
have also been reported, for example, in Refs. [ 3 , [2:]. 

In the present work, we concentrate on the produc¬ 
tion of S = —2 2 and, in particular, on the production 
reaction process of the ground state 2, 

K ( q ) + N(p) K(q') + 2 (p) , (1) 

where the arguments indicate the corresponding parti¬ 
cle’s on-shell four-momentum. This reaction has been 
studied experimentally mainly throughout the 60’s [28l — 
[37j], which was followe d by several measurements made 
in the 70’s and 80’s [38l - [44{ . The existing data are rather 
limited and suffer from large uncertainties. The total 
cross section and some of the differential cross section 
data are tabulated in Ref. [20j. We shall return to the 
discussion of these experimental data later on. Early the¬ 
oretical attempts to understand the above reaction are 
very few and can be found in Refs. [45 j-|4 9i . Recent cal¬ 
culations are reported by Sharov et al. [501 ] and by Shyam 
et al. (5l|. The former authors have considered both the 
total and differential cross sections as well as the recoil 
polarization data in their analysis, while the latter au¬ 
thors have considered only the total cross section data, 
although they too have predicted the differential cross 
sections, mentioning that they found it difficult to use 
the differential cross section data 37f for several reasons. 
Although the analyses of Refs. [50,l5l| are both based on 
very similar effective Lagrangian approaches, the num¬ 
ber of S’ = —1 hyperon resonances included in the in¬ 
termediate state are different. While in Ref. [50| only 
the E(1385) and A(1520) are considered in addition to 
the above-threshold E(2030) and E(2250) resonances^ 


2 The production threshold energy for the reaction of Eq. 0 is 


in Ref. (5l| eight of the 3- and 4-star A and E reso¬ 
nances with masses up to 2.0 GeV have been considered. 
While the authors of Ref. [Ho| pointed out the signifi¬ 
cance of the above-threshold resonances, the authors of 
Ref. M have found the dominance of the sub-threshold 
A(1520) resonance. Reaction (JT]) has been also consid¬ 
ered quite recently by Magas et al. [52j within the cou¬ 
pled channels Unitarized Chiral Perturbation approach 
in connection to the issue of determining the parameters 
of the next-to-leading-order interactions. The authors of 
Ref. [52j have added the E(2030) and E(2250) resonances 
into their calculation to improve the fit quality to the to¬ 
tal cross section data. Just recently, the Argonne-Osaka 
group [53] reported applying their Dynamical Coupled 
Channels approach to A'-induced two-body reactions for 
center-of-momentum (c.m.) energies up to IT = 2.1 GeV. 
In the reported work, both the total and differential cross 
sections were calculated, but the extracted resonance pa¬ 
rameter values are not yet available. 

We note here that the proper identification of reso¬ 
nances and the reliable extraction of their parameters 
requires detailed knowledge of the analytic structures of 
the scattering amplitude that, to date, can only be ob¬ 
tained through a full coupled-channel treatment, such as 
that of Ref. [53]. However, because the currently avail¬ 
able data in the AT2 channel are scarce and of low quality, 
they do not provide sufficient constraints for the model 
parameters to permit an in-depth analysis of that chan¬ 
nel [53j] ■ In this context, we mention that a coupled- 
channel partial-wave analysis of Af-induced reactions up 
to W = 2.1 GeV has also been performed recently by 
the Kent State University group [56[, |57[ which includes 
the KN , 7tA, ttE, ttA(1520), ttE( 1385), AT*TV, and AT A 
channels, but not the AT2 channel. 

Some of the model-independent aspects of the reac¬ 
tion ([D have been studied recently by the present au¬ 
thors [541 |55| . In the present work, we perform a model- 
dependent analysis of the existing data between thresh¬ 
old and a c.m. energy of 3 GeV based on an effective 
Lagrangian approach that includes a phenomenological 
contact amplitude which accounts for the rescattering 
contributions and/or unknown (short-range) dynamics 
that have not been included explicitly into the model. 
While the tree-level model presented here is not very so¬ 
phisticated, it captures the essential aspects of the pro¬ 
cess in question. As such, the use of a simplified, yet 
efficient model is particularly well suited for a situation, 
such as for the reaction CD, where scarce and poor data 
prevent a more detailed and complete treatment. The 
present study is our first step toward building a more 
complete reaction model capable of reliably extracting 
the properties of hyperons from the forthcoming experi¬ 
mental data, in addition to providing some guidance for 
planning future experiments. One of the purposes of the 
present work is to search for a clearer evidence of the 
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S = — 1 hyperon resonances in reaction ©■ However, we 
emphasize that our main interest here lies not so much 
in the accurate extraction of S = — 1 hyperon resonance 
parameters, but in an exploratory study to learn about 
the pertinent reaction mechanisms and, in particular, to 
identify the resonances that come out to be most relevant 
for the description of the existing S production data. In 
fact, with the exception of the £(2250) resonance, whose 
mass was adjusted slightly to better reproduce the ob¬ 
served bump structure in the total cross section in the 
charged H production, the masses and widths of the res¬ 
onances incorporated here are taken from other sources, 
as explained in Sec. ed below. Only the product of the 
coupling constants and the cutoff parameters in the cor¬ 
responding form factors are adjusted in the present work. 

The investigation of reaction © also impacts the 
study of S hypernuclei, where the elementary process of 
Eq. © is an input for the models of hypernuclei produc¬ 
tions I49l. I58M60I ] . As mentioned before, there is a pro¬ 
posed program at J-PARC and eventually at GSI-FAIR 
to obtain information about the spectroscopy of S hyper¬ 
nuclei through the antikaon-induced reactions on nuclear 
targets. Establishing the existence and properties of E 
hypernuclei is of considerable importance for a number 
of reasons and the study of reaction © is an essential 
step to this end. 

The present paper is organized as follows. In Sec. El 
our model for describing reaction © is presented, with 
some technical details supplied in the Appendix. In 
Sec. ED the results of our model calculations are pre¬ 
sented and discussed. Section Hv] contains our summary 
and conclusions. 


II. MODEL DESCRIPTION 

The reaction amplitude, T, describing a two-body pro¬ 
cess like the reaction © is, in general, given by the 
Bethe-Salpeter equation, 

T = V + VG 0 T , (2) 

where V stands for the (two-body) meson-baryon ir¬ 
reducible (Hermitian) driving amplitude and Go de¬ 
scribes free relative meson-baryon motion. Note that 
the above equation represents, in principle, a coupled- 
channels equation in meson-baryon channel space. It can 
be recast into the pole and the non-pole parts as 



T = T P + T np , 

(3) 

where the 

non-pole part T NP obeys 



t np = yNP + yN p Gq t np 

(4) 

with 




V NP = y_ V P 

(5) 


K{q) K{q') 




FIG. 1. Diagrams describing the amplitude GD in the present 
calculation. The labeling of the external legs of the s-channel 
diagram, M s , follows the reaction equation ©; the labels ap¬ 
ply correspondingly also to the external legs of the u-channel 
diagram, M u , and the contact term M c . The intermediate 
hyperon exchanges, A and £, indicated for M s also appear in 
M u - The details of the contact amplitude, M c , are discussed 
in Sec. [Ill 

denoting the one-baryon irreducible (non-pole) part of 
the driving amplitude, V. Here, V p stands for the one- 
baryon reducible (pole) part of V in the form of| 

E p = ]T|Eo r )Sor<Eo r |, (6) 

r 

where |E 0t .) and So r = (p% — m^ r + *0) -1 stand for the 
so-called bare vertex and bare baryon propagator, respec¬ 
tively. The summation runs over the baryons in the in¬ 
termediate state, each specified by the index r. The four- 
momentum and the bare mass of the propagating baryon 
are denoted by p r and m 0r , respectively. As can be seen 
in Fig. m is the sum of the s-channel Feynman dia¬ 
grams corresponding to bare baryon propagations in the 
intermediate state. The pole part of the reaction ampli¬ 
tude T p in Eq. © is given by 

T P =Y / \ F r')Sr'r{F r \, ( 7 ) 

r'r 

where the so-called dressed vertex reads 
\F r ,) = (l + T NP G 0 )\F 0r ,), 

(F r \ = (F 0r \(l + G 0 T NP ), (8) 

and the dressed propagator S r ’ r is written as 


with 

Zr'r = (F 0 r'\Go\F r ) (10) 


3 The bra and ket notation here is used only as a quick visual cue 
to identify incoming and outgoing vertices, respectively. They 
are not to be taken as Hilbert space states in the usual sense. 



4 


denoting the self-energy. 

In the present work we shall make the following ap¬ 
proximations to the reaction amplitude in Eq. ©• First, 
we approximate the pole part of the reaction amplitude 
T p by the s-channel Feynman amplitude, M s , speci¬ 
fied by effective Lagrangians and phenomenological Feyn¬ 
man propagators. Here, the dressed resonance coupling 
constants, dressed masses as well as the correspond¬ 
ing widths are parameters either fixed from independent 
sources or adjusted to reproduce the experimental data. 
The meson-baryon-baryon vertices are obtained from the 
effective Lagrangians given in the Appendix; the phe¬ 
nomenological Feynman propagators for dressed baryons 
are also found there. Note that, here, the resonance cou¬ 
plings in the dressed propagators are ignored. 

Second, the non-pole part of the reaction amplitude 
T np is approximated as follows. 

(i) Since there is no meson-exchange f-channel pro¬ 
cess in the present reaction, unless the exchanged 
meson is an exotic one with strangeness quantum 
number S = 2, f/ NP of the reaction is approxi¬ 
mated by the u-channel Feynman amplitude, M Ul 
constructed from the same effective Lagrangians 
and Feynman propagators used to construct the s- 
channel Feynman amplitudes. 

(ii) The rescattering term V NP Go T NP in T NP of 
Eq. © and other effects not explicitly included in 
the present approach are accounted for by a phe¬ 
nomenological contact term, M c , which is specified 
below. This contact term will be discussed in more 
detail later. 

With the approximations described above, the reaction 
amplitude in the present work is given by 

T = M s + M u + M c , (11) 

where M s and M u are the amplitudes from the s- and 
u-channel Feynman diagrams, respectively; both ampli¬ 
tudes include the ground-state hyperons as well as some 
of the S = — 1 hyperon resonances in the intermediate 
state. Figure |T] shows a diagrammatic representation of 
M s , M u , and M c . 

The amplitude for the reaction of Eq. 0 can be de¬ 
composed into spin-non-flip and spin-flip contributions. 
Their respective partial-wave decomposed forms reacQ 

M ++ = M__ 

= ^ E [( L + !) M l J+ (P',P) + L M T l J ~ (p',p) 

n L,T 


4 There are in total four spin matrix elements to describe the re¬ 
action 0. However, only two of them, corresponding to the 
spin-non-flip and spin-flip processes, are independent due to the 
reflection symmetry about the reaction plane for parity conserv¬ 
ing processes. See Ref. [55| for more detailed discussions. 


X P L (p-p')P T , (12a) 

M_|__ = -M_+ 

= ^ ~ Ml J -(p',p) 

n L,T 

xPl(p-p')P T , (12b) 

where initial and final momenta are as in Fig. |T] The 
indices s',s = ± in M s > s stand for spin-up (+) and spin- 
down (—) of the final (s') and initial (s) states quan¬ 
tized along the incoming momentum direction p, and 
J± = L ± \ (for L = 0, the corresponding J_ terms are 
zero). M^ J± are diagonal elements of the more general 
partial-wave amplitudes introduced in Ref. (55| (where 
full technical details can be found). The Legendre and 
associated Legendre functions are denoted by Pl(x) and 
P^ (x), respectively!! with argument p p' = cos 9, where 6 
is the scattering angle. The total angular momentum, or¬ 
bital angular momentum, and total isospin of the meson- 
baryon state are represented by J, L , and T, respectively. 
P T stands for the isospin projection operator onto the 
total isospin 0 or 1 as T = 0 or T = 1, respectively. Ex¬ 
plicitly, P T=0 = (3 + r r r 2 )/4 and P T=1 = (1 - t x ■ r 2 )/4. 

The phenomenological contact amplitude, M c , is now 
decomposed in terms of spin amplitudes similar to 
Eqs. 101) as well. Following the essential idea of Ref. j|d| , 
the corresponding contact term contributions are param¬ 
eterized as 

M c++ = M c __ 

= exp(-“z p l(p ■ P')P T , 


= J292 T (l) exp (%T Pl(p ■ p')P T , 

(13b) 

with g^ T = al exp (i^ L ), g% T = 6jexp(^ L ), and 
/3 P being constants to be fitted. (f^ L for a : = a,b 
is the complex phase angle parameter which renders the 
contract amplitude, M C1 complex and Ag is a typical scale 
parameter of the reaction at hand. The momentum de¬ 
pendence of the partial-wave matrix elements given above 
is particularly well suited for hard processes, which have 
a large momentum transfer and whose amplitudes are ex¬ 
pected to be independent of energy and nearly constant 
apart from the centrifugal barrier effects. Though reac¬ 
tion 0 is not a very hard process!! the p ,z, -dependence 


5 Here, the phase convention for the associated Legendre function 
is such that Pl(x) = sinfr). 

6 For example, the momentum transfer of this reaction at threshold 
is about 200 MeV. 
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nonetheless captures the essence of the behavior of the 
amplitude at low momentum in the final state. For fur¬ 
ther details, we refer to Ref. 0]. The exponential factor 
in Eq. (fT51) is simply a damping factor to suppress the 
high momentum behavior introduced by p' L . 

It should be noted that our phenomenological contact 
term, M c , can only account for effects with a smooth en¬ 
ergy dependence. Effects from, for example, dynami cally 
generated resonances and/or channel couplings [63465]. 
etc., that exhibit a strong variation of the amplitude as 
a function of energy cannot be described by the contact 
term. 

The amplitudes M s , M u , and M c must be added up 
to obtain the total scattering amplitude of Eq. dill) . To 
ensure that there is no ambiguity in the relative phase of 
M s + M u and M c caused by different Feynman rules, we 
give an explicit calculation of M s and M u for A(1116) in 
the Appendix. 

Standard effective Lagrangian approaches include tree- 
level s-, u- and f-channel diagrams, without phenomeno¬ 
logical contact terms. Apart from crossing symmetry de¬ 
mands, the inclusion of the u-channel amplitude, M u , in 
particular, is necessary to reproduce the backward peak¬ 
ing of the differential cross sections. (See Sec. IIIII ) In 
fact, there are a number of A and E resonances (cf. Ta¬ 
ble HI) that may contribute to this reaction. However, it 
happens that the u-channel resonance contributions, es¬ 
pecially from many of the sub-threshold resonances, also 
give rise to a total cross section which keeps increasing 
with energy in the present reaction process. This feature 
is not supported by the data, which reach a peak and 
then fall off as a function of energy. Thus, one needs a 
dynamical mechanism to suppress this rise in energy. 

The effective Lagrangian approaches of Refs. [50|, l5ll | 
have introduced phenomenological mechanisms for deal¬ 
ing with this problem that are very similar in spirit al¬ 
beit somewhat different in technical detail. In both ap¬ 
proaches, the rise of the u-channel resonance diagrams 
was suppressed with functions that smoothly cut off their 
contributions at high energies^ While the respective pro¬ 
cedures generally provide satisfactory agreement with the 
data, they both violate crossing symmetry even at the 
tree-level. 

In our model calculations, we also see the same undesir¬ 
able rise of u-channel contributions if we leave out contact 
terms. We interpret this to mean that the rescattering 
term U NP GoT NP of the non-pole T-matrix in Eq. (|4|) 
would be responsible for providing the cancellation for 
the increasing u-channel resonance amplitudes. We ac¬ 
count here phenomenologically for these in detail very 
complex dynamics by introducing contact terms, and our 
results in Sec. IIIII will show that this will indeed allow us 


7 According to a private communication by one of the authors of 
Ref. f5ll) . the form factor given in Eq. (3) of that work only 
applies to the s-channel; the u-channel was suppressed instead 
by the form factor given in Eq. (5) of Ref. [66 1. 


to treat both s- and u-channel contributions consistently, 
and at the same time avoid the high-energy u-channel 
contributions. 

In general, it seems that the problem has two scales, 
corresponding to long-range and short-range dynamics. 
The latter is, of course, sensitive to the form factors used 
at the meson-baryon vertices to account for the compos¬ 
ite nature of the hadrons, and the use of phenomeno¬ 
logical contact terms seems to be warranted to account 
for additional structure effects. Problems with two scales 
have been addressed in the past, where some authors have 
introduced two form factors, one soft and other hard, to 
mimic such effects [67]] • Also, in effective field theories 
the unknown short-range dynamics is accounted for by 
contact terms. 


III. RESULTS 

In this section, we present our results for the reaction 
K+N -A- K+E in different isospin channels. More specif¬ 
ically, we investigate the reactions K~ + p —> K + + S“, 
K~ +p —> K° +H°, and K~ + n —> K° + S - considering 
all the available data on the total and differential cross 
sections as well as recoil polarization asymmetries. 

Before we present our results, we briefly remark on 
the experimental data considered in this work, i.e., total 
cross sections, differential cross sections, and recoil po¬ 
larization asymmetries. These data come from different 
sources [dtil Kill [H, H3, EoJ and are available in various 
forms. Some of them are not in the tabular (numeri¬ 
cal) form that can be readily used but are given only in 
graphical form or as parametrization in terms of the Leg¬ 
endre polynomial expansions. In Ref. [50| , Sharov et al. 
have carefully considered the data extraction from these 
papers. We have checked that the extracted data are 
consistent with those in the original papers within the 
permitted accuracy of the check. In the present work, we 
use these data, and no cross sections resulting from the 
expansion coefficients are considered here. 

As mentioned before, there are a number of 3- and 
4-star A and E resonances, including those low-mass 
sub-threshold ones that contribute, in principle, to reac¬ 
tion ©. A list of these hyperon resonances and some of 
their properties is shown in Tabled] However, apart from 
the ground state A(1116) and E(1193), the required infor¬ 
mation for most of these resonances on the resonance pa¬ 
rameters, such as the coupling strength (including their 
signs) to S and/or N , are largely unknown. Therefore, 
the strategy adopted in this work is to consider these 
parameters as fit parameters and consider the minimum 
number of resonances required to reproduce the existing 
data. In particular, we have considered only those reso¬ 
nances that give rise to a considerable contribution to the 
cross section within a physically reasonable range of the 
resonance parameter values. More specifically, during the 
fitting procedure, resonances were added one by one to 
the model and the quality of fit was checked. It should be 



6 


TABLE I. The A and E hyperons listed by the Particle Data Group E2 (PDG) as three- or four-star states. The decay 
widths and branching ratios of higher-mass resonances ( m r >1.6 GeV) are in a broad range, and the coupling constants are 
determined from their centroid values. In the present work, the masses ( m r ) and widths (P r ) of the hyperons as given in this 
table have been used, except for the E(2250) resonance. For the latter resonance, see the text. 


A states 

E states 

State 

m r (MeV) 

T r (MeV) 

Rating 

IffiVAA'I 

State 

m r (MeV) 

r r (MeV) 

Rating 

IffiVEArl 

A(1116) 1/2+ 

1115.7 




E(1193) 1/2+ 

1193 




A(1405) 1/2" 

1406 

50 



S(1385) 3/2+ 

1385 

37 



A(1520) 3/2" 

1520 

16 








A(1600) 1/2+ 

1600 

150 

*** 

4.2 

S(1660) 1/2+ 

1660 

100 

*** 

2.5 

A(1670) 1/2" 

1670 

35 


0.3 

E(1670) 3/2" 

1670 

60 


2.8 

A(1690) 3/2" 

1690 

60 


4.0 

E(1750) 1/2" 

1750 

90 

*** 

0.5 

A(1800) 1/2" 

1800 

300 

*** 

1.0 

E(1775) 5/2" 

1775 

120 



A(1810) 1/2+ 

1810 

150 

*** 

2.8 

E(1915) 5/2+ 

1915 

120 



A(1820) 5/2+ 

1820 

80 



E(1940) 3/2" 

1940 

220 

*** 

< 2.8 

A(1830) 5/2" 

1830 

95 



E(2030) 7/2+ 

2030 

180 



A(1890) 3/2+ 

1890 

100 


0.8 

E(2250) ? ? 

2250 

100 

*** 


A(2100) 7/2" 

2100 

200 








A(2110) 5/2+ 

2110 

200 

*** 







A(2350) 9/2+ 

2350 

150 

*** 








mentioned that we have also checked the influence of var¬ 
ious combinations of resonances at a time (and not just 
one by one) to the fit quality. The resonances kept in the 
presented calculation were those that increased the qual¬ 
ity of the fit by a noticeable amount with the variation in 
X 2 per data points N, namely, &X 2 /N > 0.1. An example 
of this procedure is shown in Table Ql] where the results 
of adding one more resonance to the current model, as 
specified later, is shown. We see that some of these res¬ 
onances improve the fit quality of the total cross section 
but not the other observables or even worsen the fit qual¬ 
ity slightly. We have not included these resonances into 
our model because the total cross sections suffer from 
relatively large uncertainties. 

Whenever appropriate, for each resonance considered 
in this work, the corresponding coupling constants g KYN 
and g KY ~ were constrained in such way that the sum 
of the branching ratios Py^kn + Py^>k~. not to exceed 
unity. Because, within our model, the data are sensitive 
only to the product of the coupling constants 9kyn9kybi 
setting \g KYN \ = \9kys\ f° r the P ur POse of estimating 
the individual branching ratios, and only for this pur¬ 
pose, is a simple way of keeping our coupling constant 
values within a physically acceptable range. Admittedly, 
the currently existing data are limited and suffer from 
large uncertainties, thus an accurate determination of the 
resonance parameters are not possible at this stage. For 
this, one needs to wait for new more precise data, possi¬ 
bly including more spin polarization data. In this regard, 
the multi-strangeness baryon spectroscopy program us¬ 
ing the antikaon beam at J-PARC will be of particular 
relevance. For the ground states A(1116) and E(1193), 
the corresponding coupling constants are estimated based 
on the flavor SU(3) symmetry relations [23j |. 

It should be mentioned that, in principle, the coupling 
constants g^ Y N an( i 9kys are complex quantities ow¬ 


ing to the dressing mechanism of the resonance vertex 
as given by Eq. ([ 8 ]) . In the present work, we restrict 
them to be pure real to reduce the number of free pa¬ 
rameters and for the sake of simplicity. Note that the 
complex phases of the coupling constants are not arbi¬ 
trary in that they are constrained by unitarity of the 
scattering amplitude [H}, a feature that is absent in the 
amplitude based on a tree-level approximation as men¬ 
tioned in the introduction. However, we found that fit 
quality was not improved when we ignored theoretical 
constraints and simply allowed for complex phases in the 
coupling constants. A simple and proper way of account¬ 
ing for unitarity within an effective Lagrangian approach, 
such as the present one, is being developed and will soon 
be available elsewhere [69| ■ 

The phenomenological contact amplitude M c contains 
two sets of free parameters, {g\ T , aj} and {g % T , P^}, to 
be fixed by adjusting to reproduce the experimental data, 
for a given set of {L, T} as shown in Eq. (fl3l) . In order to 
reduce the number of free parameters, we have assumed 
the parameter aj^ to be equal to P■[ and independent 
on T and L, i.e., a^ = P J = a. The scale parame¬ 
ter A s has been fixed as Ag = 1 GeV. Note that the 
phenomenological contact amplitude can and should be 
complex in principle, since it accounts for the rescatter¬ 
ing contribution (F NP GT NP ) of the non-pole T-matrix 
which is complex in general. Accordingly, the coupling 
strength parameters g[ T and g% T are complex quanti¬ 
ties. In order to reduce the number of free parameters, 
we take their phases to be independent on L and T, so 
that, (f^ L = 4> a and p^ L = ph for all sets {L,T}. Also, 
in the present calculation, we find that it suffices to con¬ 
sider partial waves up to L = 2 in the contact amplitude 
to reproduce the existing data. 

The resonances included in the present model calcu¬ 
lations and the corresponding resonance parameters are 
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TABLE II. Variation in \ 2 per data point N, S\ 2 /N, obtained when adding one more resonance to the current model (specified 
in Table. HIIll. A negative 5x 2 /N corresponds to an improvement in the result. The quantity Sxf/Ni corresponds to Sy/ 2 /N 
evaluated for a given type of observable specified by index i = <j(total cross section), da (differential cross section), and P(recoil 
asymmetry). N = N a + N da + Np denotes the total number of data points. Furthermore, Sy 2 /Ni is given for the charged H” 
(SxL/N-) and neutral E° (Syo/No) production processes, separately. The last column corresponds to Sx 2 /N of the global fit 
considering all the data of both reaction processes. The last row corresponds to xf /Ni of the current model. 



K~ 

+ p —»• + 

£ 


I\~ 

+ p-> I< u + 




Y added 

Sxi/N a 

5XdJN da 

Sxp/Np 

<5 X-/N- 

Sxi/Nc- 

&Xd*/N da 

Sxp/Np 

Sxo/No 

Sx 2 /N 

A(1405) 

-0.01 

0.03 

0.00 

-0.01 

0.03 

0.00 

0.02 

0.01 

0.00 

A(1600) 

-0.02 

0.00 

-0.01 

-0.01 

0.02 

0.00 

0.02 

0.01 

0.00 

A(1670) 

-0.01 

0.00 

0.00 

0.00 

0.02 

0.00 

0.02 

0.01 

0.00 

A(1800) 

0.00 

0.01 

0.00 

0.00 

-0.01 

0.00 

0.01 

0.00 

0.00 

A(1810) 

-0.01 

-0.01 

0.00 

-0.01 

0.02 

0.00 

0.02 

0.01 

0.00 

A(1520) 

-0.06 

0.02 

0.00 

0.00 

-0.05 

-0.01 

0.00 

-0.02 

0.00 

A(1690) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

A(1820) 

-0.08 

0.01 

0.01 

0.00 

-0.07 

0.00 

-0.02 

-0.02 

-0.01 

A(1830) 

-0.05 

0.01 

0.01 

0.00 

0.00 

0.02 

0.02 

0.01 

0.00 

A(2110) 

-0.02 

0.02 

0.01 

0.01 

-0.03 

-0.01 

-0.03 

-0.02 

0.00 

A(2100) 

-0.08 

0.04 

0.03 

0.02 

-0.04 

-0.02 

-0.01 

-0.03 

0.01 

£(1660) 

-0.02 

0.00 

0.00 

0.00 

-0.01 

0.01 

0.00 

0.01 

0.00 

£(1750) 

-0.01 

0.01 

0.00 

0.00 

-0.01 

0.01 

0.00 

0.00 

0.00 

£(1670) 

-0.01 

0.00 

-0.01 

0.00 

0.02 

0.01 

0.01 

0.01 

0.00 

£(1940) 

0.02 

0.00 

0.01 

0.00 

0.01 

-0.01 

0.01 

-0.01 

0.00 

£(1775) 

-0.01 

0.01 

0.04 

0.01 

-0.02 

0.00 

-0.02 

-0.01 

0.00 

£(1915) 

0.01 

-0.01 

0.00 

-0.01 

-0.03 

0.00 

0.00 

-0.01 

-0.01 


Xi/Nr 

Xda/Nda 

Xp/Np 

X-/N- 

Xt/Nr 

Xda /Ndcr 

Xp/Np 

Xo/N 0 

x7 N 


1.53 

1.64 

1.89 

1.65 

0.88 

1.06 

1.74 

1.10 

1.49 


TABLE III. Fitted parameter values of the current model. For the details of the resonance parameters, see the Appendix. For 
the contact amplitude, see Eq. m- The entries in boldface are taken from Ref. [23|] and they are not fit parameters. Here, it 
is assumed that cj> PL = 4>a and (f> P = 4>b, in addition to a P = = a. 


Y 

9nak 

^NAK 

9bak 

^BAK 

A (MeV) 


A(1116)i + 

-13.24 

1.0 

3.52 

1.0 

900 








£(1193)i + 

3.58 

1.0 

-13.26 

1.0 

900 









9nak9sak 




A (MeV) 

L 

a u L 

ai 

bl 

bi 

4>a 

7 

A(1890)| + 

0.11 




900 

0 

0.28 

-1.19 





£(1385) | + 

18.76 




900 

1 

3.23 

-4.84 

-3.40 

0.61 



£(2030) l + 

0.49 




900 

2 

3.06 

21.07 

9.40 

-2.28 



£(2250)|“ 

-0.033 




900 


As = 

1 GeV 

a 

= 3.60 

0.22 

-0.16 


displayed in Table IIIII as well as the parameters of the 
phenomenological contact term, M c . We do not give 
the associated uncertainties here because they are not 
well constrained. In the present calculation, resonances 
with J < 7/2 were considered. The masses and the total 
widths of the resonances are taken to be those quoted in 
PDG na and are given in Table [TJ except for the mass of 
the £(2250) resonance. Currently, the £(2250) resonance 
is not well established and has a three-star status [la¬ 
in fact, the PDG does not even assign the spin-parity 
quantum numbers for this resonance. The analyses of 
Ref. [39^ provide two possible parameter sets, one with 
J p = 5/2~ at about 2270 ± 50 MeV and another with 
J p = 9/2“ at about 2210 ±30 MeV. In the present work 
we have assumed £(2250) to have J p = 5/2“ with the 
mass of 2265 MeV, the primary reason being that the 
total cross section in K~ + p — > K + ± S _ shows a small 


bump structure at around 2300 MeV, which is well repro¬ 
duced in our model with these parameter values. For the 
corresponding width, we have adopted the value quoted 
in PDG as shown in Table [TJ 

All parameters of the present model calculation are 
determined as described above and we now present the 
results obtained from our model. The overall fit qual¬ 
ity is quite good with x 2 /d.o.f. = 1.55 or y 2 /N = 1.49, 
as displayed in Table EH There, we also show the par¬ 
tial %-squared values xf /N i evaluated for a given type 
of observable specified by the index i as explained in the 
caption of Table EH In Fig. [2] we show the results for the 
total cross section in the charged E production reaction 
from the proton target, K~ + p —> K~ + S _ , for the 
c.m. energies up to W = 3 GeV. Figure Efa) displays 
the total contribution, which reproduces the data rather 
well. The dynamical content of the present model is also 






















FIG. 2. (Color online) Total cross section for the K~ + p —» K + + H - reaction, (a) The solid blue line represents the result 
of the full calculation of the present model. The red dashed line shows the combined A hyperons contribution. The magenta 
dash-dotted line shows the combined E hyperons contribution. The brown dotted line shows the combined A and E hyperons 
contribution. The green dash-dash-dotted line corresponds to the contact term, (b) The solid red line represents the combined 
A hyperons contribution that is the same as the red dashed line in (a). The dotted and dashed lines show the A(1116) and 
A(1890) contributions, respectively, (c) The solid magenta line represents the combined E hyperons contribution that is the 
same as the magenta dash-dotted line in (a). The dotted, dashed, dot-dashed, and dot-dot-dashed lines show the contributions 
from the E(1193), E(1385), E(2250), and E(2030), respectively. The experimental data (black circles) are the digitized version 
as quoted in Ref. j5cj] from the original work of Refs. [3614391 . 


shown in the same figure. We find that the contact term 
rises quickly from threshold peaking at around 2.1 GeV 
and falls off slowly as energy increases. It dominates the 
cross section except for energies very close to threshold 
and above ~ 2.7 GeV, where the hyperon resonance con¬ 
tributions are comparable. The A hyperons contribution 
is strongest near threshold and falls off very slowly as en¬ 
ergy increases. The E hyperons contribution is relatively 
small over the entire energy range considered, except in 
the interval of 2.0-2.3 GeV, where it becomes comparable 
to the A contribution. Near threshold, there is a strong 
destructive interference between the contact term and 
(mainly) the A hyperons contribution. At higher ener¬ 
gies, the data indicates an existence of a bump structure 
at W ~ 2.3 GeV. Our model reproduces this feature via 
delicate destructive and constructive interferences of the 
contact term and the hyperon resonance contributions. 
We also mention that we have explored the possibility 
of a much smaller contact amplitude contribution than 
shown in Fig. [2ja) considering various different sets of 
hyperon resonances from Table U however, we were un¬ 
able to find a solution with a fit quality comparable to 
that of Fig. E^a). 

Figure EKb) displays the individual A hyperon contri¬ 
butions. We see that the ground state A(1116) is, by far, 
the dominant contribution which is due to the tail of the 
corresponding u-channel process. Analogously, the indi¬ 
vidual E hyperon contributions are shown in Fig. He). 
Here, the relatively small cross section near threshold is 
due to the destructive interference between the E(1192) 
and E(1385). The enhancement of the cross section in 


the energy interval of 2.0-2.3 GeV is mostly due to the 
E(2030) resonance. The E(2250) leads to a little shoul¬ 
der in the total E contribution. We note that any non- 
negligible contribution from the hyperons for energies 
above ~ 2.3 GeV is due to the u-channel processes. 

In Fig. 02 we show the total cross section results for 
the neutral S production process, K~ + p —> K° + H°. 
Here, the data are of such poor quality that they impose 
much less constraint on the model parameters than the 
corresponding data in the charged S - production. The 
resulting dynamical content shown in Fig. 03a) is similar 
to that for the charged S - production discussed above, 
i.e., it is largely dominated by the contact term. How¬ 
ever, we see a quite different feature in the A and E reso¬ 
nance contributions as compared to that for the charged 
EG production [cf. Fig. Ha)]- One notable difference 
between the charged and neutral H production reactions 
considered here is that the u-channel A hyperon contri¬ 
bution is absent in the S° production case. Also, the 
relative contribution of the E hyperons is much larger in 
the neutral 5° production than in the charged pro¬ 
duction, especially, in the near threshold region. 

FiguresjHfb) andH(c) show the individual hyperon con¬ 
tributions. As mentioned before, due to the absence of 
the u-channel A exchange in the neutral £° production, 
the A(1116) contribution is insignificant, leading to a neg¬ 
ligible contribution of the A hyperons. Due to the isospin 
factors here, the E(1192) and E(1385) hyperons interfere 
constructively, especially near the threshold. Recall that, 
for charged 5 - production, these hyperons interfere de¬ 
structively [cf. Fig. Hr)]- 
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FIG. 3. (Color online) Same as Fig. [5] for the K + p —» K° 
digitized version as quoted in Ref. [50|] from the original work c 



+ S° reaction. The experimental data (black circles) are the 
f Refs. 




FIG. 4. (Color online) Total cross section results with individual resonances switched off (a) for K~ + p —» K + + S“ and (b) 
for K~ + p —> K° + E°. The blue lines represent the full result shown in Figs. [2] and QT] The red dashed lines, which almost 
coincide with the blue lines represent the result with A(1890) switched off. The green dash-dotted lines represent the result 
with E(2030) switched off and the magenta dash-dash-dotted lines represent the result with E(2250)5/2“ switched off. 


In Fig. H] we illustrate the amount of the above¬ 
threshold resonance contributions of the present model 
to the total cross sections. We do this by comparing the 
full results shown in Figs. 0(a) and O a ) to the result 
found by switching off one resonance at a time. We see 
in Fig. 0K a) that the largest effect of £(2030) on the cross 
sections is in the range of W ~ 2.0 to 2.4 GeV. This res¬ 
onance is clearly needed in our model to reproduce the 
data. It also affects the recoil polarization as will be dis¬ 
cussed later. It should be mentioned that this resonance 
also helps to reproduce the measured K + E~ invariant 
mass distribution in 7 + p —»• K + + K + + 5“ (24|, by 
Filing in the valley between the two bumps in the in¬ 
variant mass distribution that would appear without it; 
such a feature clearly is not observed in the data 0 . 
The A(1890) affects the total cross section in the range 
of W ~ 1.9 to 2.1 GeV, and the £(2250)5/2“ contributes 


around W ~ 2.2 GeV, where it is needed to reproduce 
the observed bump structure. A more accurate data set 
is clearly needed for a more definitive answer about the 
roles of the A(1890) and £(2250) resonances. FigureHKb), 
for the neutral 5° production, also shows a similar fea¬ 
ture observed in the E“ case for the £(2030) resonance. 
Here, the influence of the £(2250)5/2“ is smaller and 
that of the A(1890) is hardly seen. Recall that there is 
no u-channel A contribution in the neutral 5° produc¬ 
tion. 

The results for differential cross sections in both K~ + 
p —> K + + E“ and K~ + p —> K° + S° are shown in 
FigS-EKa) and(5](b), respectively, in the energy domain up 
to W = 2.8 GeV for the former and up to W = 2.5 GeV 
for the latter reaction. Overall, the model reproduces the 
data quite well. As in the total cross sections, the data 
for the neutral 5° production are fewer and less accurate 
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K+p^K + +E (a) K+p^K°+E° (b) 




FIG. 5. (Color online) Kaon angular distributions in the center-of-mass frame (a) for K~ + p — » K + + E _ and (b) for 
K~ + p — > K° + S°. The blue lines represent the full model results. The red dashed lines show the combined A hyperons 
contribution. The magenta dash-dotted lines show the combined E hyperons contribution. The green dash-dash-dotted line 
corresponds to the contact term. The numbers in the upper right corners correspond to the centroid total energy of the system 
W. Note the different scales used. The experimental data (black circles) are the digitized version as quoted in Ref. |50| from the 
original work of Refs. [3ll-f34l. l36i. [37] for the K~ + p — > K + +E _ reaction and of Ref. "30, 36, 43, ’43 for the K~ +p —» K° + E° 
reaction. 



S wave 


W (GeV) 


K +p^ K + + S (b) 



FIG. 6. (Color online) Partial wave decomposition of the total cross section and the angular distribution for K~ +p —> K + +E“. 
(a) Total cross section sectioned by contributions from each partial wave L. The red shaded area indicates the S -wave 
contribution, while the green area corresponds to the P-wave. Magenta indicates the D -wave and maroon the F- wave, (b) K + 
angular distribution: the solid blue lines are the full results, while the dotted green lines represent the sum of S + P waves, the 
red dashed lines represent the S + P + D waves and the dash-dotted magenta lines represent the S + P + D + F waves. For 
lower energies, the S + P + D waves already saturate the full cross section results so that the F- and higher-wave contributions 
cannot be distinguished from the full result. 


than for the charged 2“ production. In particular, the 2° 
production data at W = 2.15 GeV seems incompatible 
with those at nearby energies, and the present model is 
unable to reproduce the observed shape at backward an¬ 
gles. It is clear from Figs.Oa) and [51(b) that the charged 


channel shows a backward peaked angular distributions, 
while the neutral channel shows enhancement for both 
backward and forward scattering angles (more symmet¬ 
ric around cos 0 = 0) for all but perhaps the highest ener¬ 
gies. In the charged 2 _ production, both the resonance 
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FIG. 7. (Color online) Ratio of the measured total cross sec¬ 
tion a and the final state A'E relative momentum p' as a func¬ 
tion of p' 2 . The blue square data correspond to K~ + p —» 
A + + E“, while the red circle data to K~ +p —¥ A^ + E 0 . The 
blue solid and red dashed curves are the present model results 
corresponding to K~ +p —> K + + E _ and K~ +p —» K° + E°, 
respectively. 


and contact amplitude contributions are backward angle 
peaked and, as the energy increases, they get smaller and 
smaller at forward angles. In 2° production, both the £ 
resonance and contact amplitude contributions also ex¬ 
hibit an enhancement for forward angles. Note that the 
A resonance contribution here is negligible due to the 
absence of the u-channel process. The interference pat¬ 
tern in the forward angular region depends on energy. 
At lower energies the interference is constructive and it 
becomes destructive at higher energies. The behavior of 
the angular distributions in terms of the partial waves 
will be discussed later in connection with the results of 
Figs. [(2b) and [Hb). 

The partial-wave content of the cross sections for the 
charged 2“ production process arising from the present 
model is shown in Figs. Eta) and (Htb). As can be seen 
in Fig. Ha), the total cross section is dominated by the 
P and D waves in almost the entire energy range con¬ 
sidered, even at energies very close to threshold where 
one sees a strongly rising P -wave contribution. The S- 
wave contribution is very small. This peculiar feature is 
caused by the ground state A(1116), whose contribution 
cancels to a large extent the otherwise dominant S'-wave 
contribution close to threshold, in addition to enhancing 
the P -wave contribution. One way of probing the S-wave 
content close to threshold in a mo del-independent man¬ 
ner is to look at the quantity a/p' as a function of p' 2 , 
where p' is the relative momentum of the final A'E state. 
The reason being that, for hard processes, the partial- 
wave reaction amplitude behaves basically like p' L for 
a given orbital angular momentum L as mentioned in 
Sec. [HI This leads to 

^7 = c o + Ci p' 2 + c 2 p' 4 + ... , (14) 


with expansion constants c L . Figure [7] illustrates this 
point. Although the existing experimental data are of 
poor quality, they reveal the general features just men¬ 
tioned. In particular, for the charged 2“ production pro¬ 
cess, the data indicate a linear behavior of a/p’ close to 
threshold implying a strong P -wave contribution. Our 
present model result is consistent with this behavior. It 
is also consistent with the observation made in Ref. (53j 
that the low-energy behavior of the total cross sections 
in the ttN, r]N and KE channels does not seem to follow 
the usual S'-wave dominance. 

The corresponding results for the neutral 2° produc¬ 
tion are also shown in Fig. [7] There, the scattered 
data are consistent with S-wave dominance, a feature 
exhibited by our model as well [see also Fig. [8) a )]- hi 
Fig. [6D a )i we also see a small F -wave contribution above 
W ~ 2.0 GeV that helps saturate the total cross section. 
Note that since our contact term includes partial waves 
only up to L < 2, the F -wave contribution is entirely 
due to the hyperon resonances. The enhancement of the 
D -wave contribution around W = 2.3 GeV as well as 
the little shoulder in the P -wave contribution are due to 
the £(2250) hyperon. Of course, the partial-wave contri¬ 
butions are mainly constrained by the differential cross 
section and they are shown in Fig. [Gjb) . As mentioned 
before, the shape of the angular distribution is backward- 
angle peaked and the cross sections are very small at for¬ 
ward angles. This behavior is a direct consequence of the 
very significant interference between the P and D waves. 
This can be seen by expanding the cross section in partial 
waves. Considering the partial waves through L = 2 and 
following Ref. [55], the differential cross section may be 
expressed as 


da 

dfl 


*02 I 


|aq| + 2Re(a 02 d 2 ) 


cos 2 9 


+ |a 2 | cos 0 + 


8 + (l/3i| 2 + l/3 2 | 2 cos 2 sin 2 8 


2 Re 


a 02 oq + oi x a 2 cos 2 9 + /3 x jd 2 sin 2 6 


cos 8 


(15) 


where the coefficients a L (/3 L ) provide a linear combina¬ 
tion of the partial-wave matrix elements corresponding 
to the spin-non-flip (spin-flip) process with a given or¬ 
bital angular momentum L [see Eq. (1121) ]. Here, a 02 = 
a 0 — |d 2 , a 2 = |a 2 , and /3 2 = 3/3 2 . In the above equa¬ 
tion, the last term on the right-hand side involving an 
interference between the P and D waves is an odd func¬ 
tion in cos 0, while the the first term in square brackets 
is an even function. These two terms cancel to a large 
extent at forward angles while at backward angles they 
add up. Note that these partial waves are comparable in 
strength as shown in Fig. EJa) so that their interference 
term leading to an odd function part can largely cancel 
the even term at forward angles. 

Figures and (Hlb) display the partial wave con¬ 
tent in the cross sections for the neutral 2° production 
process. In contrast to the charged 2“ production, here 
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K + p —> K°+ Z° (b) 




FIG. 8. (Color online) Same as in Fig. ©but for K~ + p -» K° + S°. 


one sees that the largest contribution to the total cross 
section is the D- wave, and the P- wave is largely sup¬ 
pressed, which is a direct consequence of the shape of the 
observed angular distribution whose partial wave contri¬ 
butions are shown in Fig. [8{b). There, compared to that 
for charged 2“, one sees a more symmetric angular shape 
about cos 9 = 0 that is dominated by the D- wave. The 
present model reproduces the observed behavior of the 
K° angular distribution by suppressing the P- wave con¬ 
tribution as can be easily understood from Eq. (1151) . The 
rather drastic suppression of the P wave can be better 
seen in Fig. [81(a). For energies very close to threshold, 
the cross section is dominated by the S- wave as seen also 
in Fig. [HI a). 

The results for the recoil polarization asymmetry mul¬ 
tiplied by the cross section are shown in Fig. [9] in the 
energy interval of W = 2.1 to 2.5 GeV. Overall, we 
reproduce the data reasonably well. We also find that 
the results shown at W = 2.11 GeV are still signifi¬ 
cantly affected by E(2030). This corroborates the find¬ 
ings of Ref. [50:] . We recall that the recoil asymmetry 
is proportional to the imaginary part of the product of 
the non-spin-flip matrix element (M ss ) with the com¬ 
plex co njug ate of the spin-flip matrix element ( M s i s with 
s' ^ s) J55(, so that it vanishes identically unless these 
matrix elements are such that their product has a non¬ 
vanishing imaginary part. We can therefore expect the 
recoil polarization to be sensitive to the complex na¬ 
ture of the reaction amplitude, in particular, to the phe¬ 
nomenological contact amplitude, M c , introduced in the 
present model. Indeed, if one forces the coupling strength 
parameters, <?f T and T in Eq. (fT3l) . to be pure real dur¬ 
ing the fitting procedure, the Xp/Np deteriorates, e.g., 
from 1.89 to 2.26 for the K~ + p —> K + + 2“ reaction, 
although the quality of fit for cross sections is nearly un¬ 
changed. 


K + p —> K + K + p —> K + cs 



cos 0 cos 9 


FIG. 9. (Color online) The recoil asymmetry multiplied by 
the cross section, P for both the K~ +p — > K + + H - and 
K~+p —¥ K°+E° reactions. The blue solid lines represent the 
full results of the current model. Data are from Refs. [dll. f37f|. 


In Fig. 1101 we show the present model predictions for 
the target-beam asymmetries, K xx and K xz , multiplied 
by the unpolarized cross section, i.e., j^K xx and ^ K xz 
for both the charged EV and neutral 2° production pro¬ 
cesses. These observables are related to the spin-rotation 
parameter /3 [70j by tan/3 = —K xz /K xx . Note that 
these target-recoil asymmetries, together with K yy , are 
the only three independent double-spin observables in the 
reaction of Eq. ([]]) as discussed in Ref. [55} . Indeed, the 
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K +p^> K + + E (a) 


K +p-> K°+ E° (b) 



FIG. 10. (Color online) Target-recoil asymmetries K xx (green dashed curves) and K xz (blue solid curves) as defined in Ref. f55t ] 
multiplied by the cross sections for the reactions (a) K~ + p —» K + + H _ and (b) K~ + p —> K° + S°. The numbers in the 
upper right corners represent the total energy of the system W in units of GeV. 


only two other non-vanishing target-recoil asymmetries 
are related by K zz = K xx and K zx = —K xz @ We men¬ 
tion that ^ K xx is proportional to the difference of the 
magnitude squared of the spin-non-flip and spin-flip ma¬ 
trix elements, while j^K xz is proportional to the real 
part of the product of the spin-non-flip matrix element 
with the complex conjugate of the spin-flip matrix ele¬ 
ment. Therefore, unlike the recoil asymmetry, these spin 
observables do not vanish even if the reaction amplitude 
is pure real or pure imaginary. This means that they are, 
like the cross section, much less sensitive to the complex 
nature of the phenomenological contact amplitude. 

To gain some insight into the angular dependence ex¬ 
hibited by these target-recoil asymmetries in Fig. [TUI we 
express them in terms of partial waves with L < 2, which 
gives 


dn Kxx ~ |a ° 21 


|aq| + 2 Re (a 02 d 2 ) 


cos 2 9 


+ |d 9 | cos 4 9 — 


II t |/3 2 | 2 COS 2 0 ) sin 2 0 
2 a 


da 

dQ 


K xz = 2 Re 


2 Re a 02 aq + aqa 2 cos^ 

—/3 1( 3 2 sin 2 0 cos0 , (16a) 

*02/3* + ^ 1 P 2 + cos 2 9 


8 Note that the symmetry of the reaction leads to K yy = 7r- inde¬ 
pendent on the scattering angle 6 [54l [HU] . Here, 7 r w stands for 
the parity of the produced S which is taken to be 7 t ct = +1 for 
the ground state S. Also, K xx = K zz | cog = 7r w . The tar¬ 
get asymmetry is identical to the recoil asymmetry in the present 
reaction. Therefore, we exhaust all the independent observables 
available in the reaction processes considered here. 


( a 02$2 + a lPl) cos 0 + U 2/ d 2 cos 3 9 


sm( 


(16b) 


Note that the only difference between j^K xx given above 
and differential cross section given by Eq. m is the sign 
change of the terms involving (3 L . These terms are, how¬ 
ever, proportional to sin 2 9. Therefore, this spin observ¬ 
able behaves like the differential cross section at very 
forward and backward angles, where sin 2 9 -C 1. At 
cos 9 = 0, the difference is due to the term of ±|/3i| 2 , 
which is a P- wave contribution in the spin-flip amplitude. 
Now, if we ignore the P- wave contribution — which is rel¬ 
atively very small in the neutral S° production over the 
nearly entire energy region considered as seen in Fig.[8la) 
it is immediate to see that Eq. (TTUh ) involves only 
terms that are symmetric about cos0 = 0. We see in 
Fig. Mb) that j^K xx exhibits roughly this symmetry. 

For j^K xz , Eq. (fT6b) reveals a rather complicated 
angular dependence in general, and no particular fea¬ 
ture is apparent in the results shown in Fig. 1101 espe¬ 
cially for the charged H - production process. Neglecting 
the P- wave contribution, Eq. (fTUb) reduces to %^K XZ = 


Re 


(a 02 + a 2 cos 2 9) f3* 


sin 29, which is roughly the an- 


02 ~ u 2 LUD u J r'2 
gular dependence exhibited in Fig. [TOl' b). 

The present model predictions for the K~ + n —> K° + 
S - reaction are shown in Fig. llll Here, the experimental 
data are extremely scarce, and they were not included in 
the present fitting procedure. Nevertheless, the current 
model is seen to predict those few data quite reasonably. 
Both the total and differential cross sections exhibit a 
very similar feature to those of the K~ + p —> K + + S - 
reaction with a noticeable small enhancement in the dif¬ 
ferential cross sections as seen in Fig. [TTT b) for forward 
angles near cos 9 = 0 in K~+n —> K°+'E~. We see, how¬ 
ever, some bigger differences in the individual amplitude 
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K + n -> A°+ S (a) 



FIG. 11. (Color online) Same as Figs. [2! a) and [5] for the K 
Refs. [HIM]. 


contributions, more clearly seen in the total cross sections 
that are given in Fig. [TlT a). There, the E hyperon con¬ 
tribution is larger than the A contribution over the entire 
energy region up to W ~ 2.3 GeV, in particular, at low 
energies near threshold. This is due to the absence of the 
strong destructive interference between the E(1385) and 
E(1192) (not shown), since the latter hyperon contribu¬ 
tion is suppressed to a large extent compared to the case 
of K~ +p —> K + + 2 _ . Moreover, there is a constructive 
interference with the A hyperon, which makes the sum 
of the hyperons contribution relatively large in the low 
energy region. 

For completeness, we also show in Fig. [12] results for 
the Kl + p —> K + + 2° reaction. Within the present 
model, the cross sections for this process simply differ 
by a factor of 1/2 from those shown in Fig. [FT] for the 
K~ + n —> K° + 2“ reaction. We show the Kl results 
here because the creation of a high-intensity Kl beam 
currently being contemplated 0 may open up an entire 
new and exciting field of hyperon spectroscopy. 


IV. CONCLUSION 

In this work we have presented our analysis on the re¬ 
action of K~ + N —> AT+2 within an effective Lagrangian 
approach that includes a phenomenological contact term 
to account for the hnal-state-interaction rescattering con¬ 
tribution of the reaction amplitude in the Bethe-Salpeter 
equation and for other possible (short-range) dynamics 
that are not explicitly taken into account in the model. 
By introducing this phenomenological contact term, we 
avoid the problems found in the usual implementations 
of tree-level effective Lagrangian approaches that need 
to phenomenologically suppress u-channel contributions 
dominating the high-energy behavior j50l l5lj] . In addition 
to the ground states A(1116) and E(1193), the present 



cos 0 


+ n —> K° + S reaction. The experimental data are from 


model also includes the A(1890), E(1385), E(2030), and 
E(2250) resonance contributions. 

The available total and differential cross sections, as 
well as the recoil asymmetry data, in both the K~ +p —> 
K + +2“ and K~ + p —> K° + 2° processes are well 
reproduced by the present model. We have found that 
the above-threshold resonances A(1890), E(2030), and 
E(2250) are required to achieve a good fit quality for the 
data. Among them, the E(2030) resonance is the most 
critical one. This resonance affects not only the cross sec¬ 
tions but also the recoil asymmetry. In addition, it also 
brings a model calculation of Ref. [24] into an agreement 
with the observed 7F + 2 _ invariant mass distribution in 
2 photoproduction [l9|. The A(1890) is also required to 
improve the fit quality in the present model, especially 
in the energy dependence of the total cross sections of 
the charged 2“ production around W = 1.9 GeV. The 
total cross section data in the charged 2 _ production 
seems to indicate a bump structure at around W = 2.3 
GeV, which is accounted for by the E(2250) resonance 
with J p = 5/2— and a mass of 2265 MeV in the present 
model. More accurate data are required before a more 
definitive answer can be provided for the role of these two 
resonances. In this regard, the multi-strangess hyperon 
production programs using an intense antikaon beam at 
J-PARC is of particular relevance in providing the much 
needed higher-precision data for the present reaction. 

The present analysis also reveals a peculiar behavior 
of the total cross section data in the threshold-energy re¬ 
gion of the K + 2“ production channel, where the higher 
partial-waves (P and D) dominate instead of the usual 
S'-wave (cf. Fig. [7]). If this behavior of the cross sec¬ 
tion data is corroborated in future experiments, it will 
cast serious doubts on the validity of model calculations 
that neglect higher partial-waves even for energies very 
close to threshold. This peculiar low-energy behavior of 
the total cross section in AT-induced reaction seems to be 
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present also in the ttN and ttN production channels, in 
addition to the A'3 channel [53[ ■ 

Apart from the recoil asymmetry, we have also pre¬ 
dicted the target-recoil asymmetries for which there are 
currently no experimental data. In contrast to the re¬ 
coil polarization — which are small these observables 
are quite sizable and may help impose more stringent 
constraints on the model parameters. In principle, one 
requires the four independent observables calculated here 
to completely determine the reaction amplitude [55| . Of 
course, measurements of the spin obervables, in particu¬ 
lar, are challenging experimentally by any standard, but 
one may exploit the self-analyzing nature of the produced 
hyperon to help extract these observables [55|, uM • F° r 
the target-recoil asymmetry measurements, one requires 
a polarized target in addition to spin measurements of the 
produced 5. Polarized targets available at some of the 
world’s major laboratories combined with the availability 
of intense beams make measuring these spin observables 
no longer out of reach. In fact, various single- and double¬ 
polarization observables in photoproduction reactions are 
currently being measured at major facilities such as JLab, 
ELSA, and MAMI, aiming at so-called complete experi¬ 
ment sets in order to model-independently determine the 
photoproduction amplitudes. 

While it may perhaps not be entirely clear which role 
any particular resonance plays for the K~ + N —> K + E 
reaction, the present and other calculations based on ef¬ 
fective Lagrangians da nn , and also the unitarized chi¬ 
ral perturbation approach [52|, seem to agree that some 
S = — 1 hyperon resonances seem to be required to re¬ 
produce the existing data. To pin down the role of a 
particular resonance among them requires more precise 
and complete data, in addition to more complete theo¬ 
retical models. In any case, the present reaction is very 
well suited for studying S = — 1 hyperon resonances. 

Finally, the present work is our first step toward build¬ 
ing a more complete reaction theory to help analyze the 


data and extract the properties of 3 resonances in fu¬ 
ture experimental efforts in 3 baryon spectroscopy. This 
is a complementary work to that of a model-inde pen dent 
analysis performed recently by the same authors [53 and 
will also help in analyzing the data to understand the 
production mechanisms of 3 baryons. 
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Appendix 

In this Appendix, we give the effective Lagrangians 
and phenomenological dressed baryon propagators from 
which the s- and u-channel amplitudes, M s and M u dis¬ 
cussed in Sec. El are constructed. We follow Refs. QEl 
Ito! [75l | and consider not only the spin-1/2 ground state A 
and E but also their respective excited states with spin 
up to 7/2. In the following we use the notations for the 
iso-doublet fields 
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and for the iso-triplet fields 


where m K denotes the kaon mass. For spin-5/2 hyperons 


£ = 



(A.2) 


We also introduce the auxiliary operators in Dirac 
space 

= -r (±) (±»A + l ~ X $) , (A.3a) 

\ rn B' ^ rn B J 

D 3/2(±) = T (T)d„ 7 (A.3b) 

D% 2(±) = , (A.3c) 

(A.3d) 


where = y 5 and Tl - ) = 1. Here, m B stands for 
the mass of the baryon B. The parameter A has been 
introduced to interpolate between the pseudovector (A = 
0) and the pseudoscalar (A = 1) couplings. Note that in 
the above equation the order of the subscript indices in 

rj 1 / 2 ^) • ; 1Tmor f an t j p jj 1 / 2 ^) / rj 1 / 2 ^) 

b'bm 1S lm porrant, i.c., r 1 j bb'm- 

The effective Lagrangians for spin-1/2 hyperons A and 


£ (or their resonances 

:) are, then, given by 


pi/ 2 (±) _ „ x( 

t-ANK ~ 9ANK l 

DankK)n + H.c., 

(A.4a) 

^>l/2(±) _ y, 

l "snk — 9y,nk ^ 

(dTnkK) tAT + H.c., 

(A.4b) 

cV 2 (±) _ _ u( 

l -eak c — 9sak c - ( 

d sak )k c) A + H-c., 

(A.4c) 

l -ey,k c — 9st,k c - T 

■ ( d svk )k c) ■ s + H.C., 

(A.4d) 

where the superscripts ± refer to the positive (+) and 
negative (—) relative parity of the baryons. Flavor 


SU(3) symmetry relates the coupling constants among 
the members of the octet J p = l/2 + ground state 
baryons and J p = 0~ pseudoscalar mesons and we have 


1 + 2a 

9ank — ~98 — yjj— ’ 

9y.nk = 38 (1 — ^ a ), 
1 — 4a 

9sak c ~9& i 

9ssk c = —9 8 i 


(A.5a) 
(A.5b) 
(A.5c) 
(A.5d) 


where the empirical values are g s = g^Nn = 13-26 and 
a = 0.365, where a is the F/D mixing parameter defined 
as a = F/(D + F). 

For spin-3/2 hyperons, we have 

£ 3/2(T) = 9^k_ ^ D 3/2 (±)^ N + H c t (A.6a) 

£ 3/2(±) = g^K_ . ^ D 3/2(±)x>j tN + H c 5 (A. 6b ) 

5 [D 3 J 2 ^K^j A" + H.c., (A.6c) 

44 2 S } = ( D l /2(±)K c)'+ H.c., (A.6d) 


4 !nk = ^F ~^ V (d%^k) N + H.c., (A.7a) 

C VNK = ^F ^ ■ (d 5 J? (±) k) tN + H.c., (A.7b) 
4 a£ } = S (d%^K c ) AT + H.c., (A. 7c) 
C %Sr ( D% 2 ^K c ^ ■ + H.c.. 

(A.7d) 

And, for spin-7/2 hyperons, we have l2~il l~76l] 

^Ink = ^F ~^ p ( d I /2 p (±) k) n + H-c., (A.8a) 

&k = ^ ■ ( d I /2 p (±) k) tN + H.c., 

(A.8b) 

4® = S (d% 2 p {±) K c ) ATP + H.c., (A.8c) 

(dICT^Kc) ■ ^ vp + H.c.. 

(A.8d) 


The coupling constants in the above Lagrangians cor¬ 
responding to A and £ resonances are free parameters 
adjusted to reproduce the existing data. For those reso¬ 
nances considered in the present work, they are given in 
Table ITITI 

In the present work, all the meson-baryon-baryon ver¬ 
tices are obtained from the above Lagrangian. In addi¬ 
tion, each vertex is multiplied by an off-shell form factor 
given by 


/ (j?r , , Ar ) — 


nAt 


nAf. + (p 2 — to)!) 2 


(A.9) 


where p 2 and m r are the square of the 4-momentum and 
mass of the exchanged hyperon, respectively. The cutoff 
parameter A r is chosen to have a common value A r = 
A = 900 MeV for all the MBr vertices in order to keep 
the number of free parameters to a minimum. Also, we 
choose n = 1. 

For the propagators of the dressed hyperons, we 
could in principle adopt the forms used in our previous 
work [24l. I73IT75| . However, in view of the limited amount 
of currently available data for the present reaction and 
the rather poor quality of these data, here we adopt the 
simpler forms as given in the following. For a spin-1/2 
baryon propagator, we use 


sy 2 (pr) 



(A.10) 


where r r is the baryon width assumed to be constant, in¬ 
dependent of energy. For a stable (ground state) baryon, 
r r e with e being positive infinitesimal. 
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For spin-3/2, the dressed propagator reads in a 
schematic matrix notation 

Sl /2 { Pr ) = - --A , (A.11) 

t r - ™ r + 

where A is the Rarita-Schwinger tensor with elements 


A yv +44 

O OTll r 


r y^ L p 1 ' — 

3 m„ 


. (A.12) 


Similarly, the propagator for a spin-5/2 resonance is given 

by 


S 5 r /2 (Pr ) = 


— m r + i- 


-A. 


(A.13) 


where the elements of A are 

. ft. 1 


a l ° c 2 2 


( ga\gil +gi\gt\} - ^g ai a 2 g ^ 2 

(ja^g'k+i'cxj^Sal +7a 2 7 /3l 5a 2 1 


1 

10 


+ 7a 2 7 ^gi\) 


(A.14) 


with 


r v = <r 


■p^p” 


(A.15) 
mi. 


The propagator for a spin-7/2 resonance is given by 

1 


S 7 r /2 (Pr) = 


— ro„ + i- 


-A. 


(A.16) 


where the elements of A are 


a “l°2«3 / ■, 


P(a),P(P) 


—y I —H2 —A'S — /3i — —/5 

ga 1 ga 2 ga 3 — jg ai g a2 a 3 g 


_ PyP 1 Qp2 I Q qP2P3 

pa.fl ga 2 ga 3 -T 3g7 ai 7 y<x 2 a 3 g 


(A.17) 

and the summation runs over all possible permutations 
of {anaj.aj} and of {A,/3 2 ,/? 3 }. 

To avoid an ambiguity in the relative phase between 
M s + M u and M c in Eq. (ED, we provide here the ex¬ 
plicit expressions for the amplitudes M s and M u for the 
A(1116) exchange in the K~{q) + N(p) —1 K + (q') + 
S ~{p r ) reaction, i.e., 

M s a = , 

(A.18a) 

M„ a = u 3 (p f ) r l- K - K {q) S]! 2 { Pu ) r U AK+N (q') u N (p ) , 

(A.18b) 

where the nucleon index N stands for the proton and 
A stands for A(1116); the baryon Dirac spinors are nor¬ 
malized covariantly, u B u B = 1, the intermediate four- 
momenta are p s = p + q and p u = p — q', and the vertices 
are given as 

t ak-n(<i) = 9ank 7 5 (a - 1 X i) fs , (A.19a) 

x m A+p i N y 

^k + aW) = 9eak 7 s A + $ f s , (A.19b) 

\ 11 I I It J 

^ak+nW) = 9ank 7 5 (a + — f u , 

V m A ”T m N / 

(A.19c) 

V Ik- a(«) = 9sak 7 5 ^A - • fu , (A.19d) 

where A describes the linear interpolation between pseu¬ 
doscalar (A = 1) and pseudovector (A = 0) couplings; the 
off-shell form factors are given by [see Eq. (IA.9I) ] 

fs = /(Ps,mA,A A ), (A.20a) 

fu = f{p 2 u, m A , Aa) , (A.20b) 

and the values of Panki 9~aki an d A are found in Ta¬ 
ble Ell 
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